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Abstract. We analyze a finite volume scheme for nonlocal SIR model, which 
is a nonlocal reaction-diffusion system modeling an epidemic disease. We es- 
tablish existence solutions to the finite volume scheme, and show that it con- 
verges to a weak solution. The convergence proof is based on deriving series 
of a priori estimates and using a general L p compactness criterion. 



1. Introduction 



We consider a mathematical model describing an epidemic disease in a physical 
domain ft C R d (d = 1,2,3) over a time span (0, T), T > 0. In this model 
we consider the propagation of an epidemic disease in a simple population p = 
u\ + U2 + Us, where u\ = u\(t, x), U2 = U2(t, x) and 1*3 = Ua(t, x) are the respective 
densities of susceptible (those who can catch the disease), infected (those who have 
the disease and can transmit it) and of recover individuals (those who have been 
exposed to the disease and will become infective after the lapse of an incubation 
period) at time t and location x. A prototype of a nonlinear system that governs 
the spreading of a nonlocal SIR of epidemics with through a population in a spatial 
domain is the following nonlocal reaction-diffusion system: 

dtU\ — <zi (^J u\ dxj Aui — —a(ui,v,2, U3) — /xui, 

(1.1) < d t U2 - Ctef / U2dx J AU2 = a(ui,U2,U3) - 7K2 - /XH2, 

in Qt, where Qt denotes the time-space cylinder (0,T) x fL We complete the 
system (jl.lj) with Neumann boundary conditions: 



(1.2) 



u, dx^j Vu, ■ T) = ondOx(G,T), i = 1,2,3, 



where 77 denotes the outer unit normal to the boundary d£l of Q, and with initial 
data: 



(1.3) 



1,2,3. 
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One of the simplest SIR models is the Kermack-McKendrick model [TT] which 
consists in a systems of 3 x 3 EDOs which is similar to the model above, but 
without diffusion terms. The classical epidemic SIR model describes the infection 
and recovery process in terms of three ordinary differential equations has been 
studied by many researchers [18] and the reference cited therein. 

Hence the description of the classical SIR model in the system (jl.ip . i is the 
length of latency period or duration of the exposed stage, and /i the natural mor- 
tality rate. The incidence term a take the following form: 

uv 

(1.4) c(w, v, w) — a for some a > 0, 

u + v + w 

which coincides with the classical model cr(u, v, w) — auv when the total population 
P(t) remains constant. In fact, the well known SIR model which appears generally 
in literature is renormalized which we do not suppose here. For technical reasons, we 
need to extend the function a(u,v,w) so that it becomes defined for all (u,v,w) £ 
Ixlxl. We do this by setting 

a(u+, v+, w+) if (u, v, w) ^ (0, 0, 0), 

otherwise . 



a(u, v, w) = 



In this work, the diffusion rates a\ > 0, a 2 > and 03 > are supposed to 
depend to the whole of each populations in the domain rather than on the local 
density, i. e. moves are guided by considering the global state of the medium. We 
assume that ai : R — > R is a continuous function satisfying: there exist constants 
Mi, C > such that 

(1.5) Mi < a. t and |a;(/i) - a t {I 2 )\ < G \h - I 2 \ for all I U I 2 £ R, for i = 1,2, 3. 

Such equations with nolocal diffusion terms has already been studied from a theo- 
retical point of view by several authors. First, in 1997, M. Chipot and B. Lovat [1] 
studied the existence and uniqueness of the solutions for a scalar parabolic equation 
with a nonlocal diffusion term. Existence-uniqueness and long time behaviour for 
other class of nonlocal nonlinear parabolic equations and systems are studied in 
[IJ[14]. Liu and Jin made some experimental simulations in 13 in order to observe 
spacial patterns in an epidemic model with constant removal rate of the infective. 

Before we define our finite volume scheme, let us state a relevant definition of a 
weak solution for the nonlocal SIR model. 

Definition 1.1. A weak solution of Ql.l[) - (|1.3[l is a triple u = (iti, U2, W3) of func- 
tions such that Mi, it 2 , U3 E £ 2 (0, T; H 1 ^)), 

(1.6) 

Uifi(x)ipi (0, x) dx — jj u\dt^\dxdt+ II a i[ u x dx I Viti • V^i dxdt 
(a(ui,u 2 ,u 3 ) + nui)ip 2 dxdt, 



(1.7) 

— / M 2. 0(^)^2(0, x) dx — 1 1 u 2 dt<p2 dx dt + // a 2\ u 2 dx]\7u 2 - Vip 2 dxdt 
Jn JJqt J JQt ^Jn ' 



(a(ui,u 2 ,u 3 ) - (7 + fi)u 2 )ip 2 dxdt, 
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(1-8) 

it3 dx ) V1/3 • Vif3 dx dt 



^3,0(^)^3(0, x) dx - // u 3 d t tp3dxdt + // a 3 ( 

JJqt J Jq t Vn 

7«,2(/?3 dxdt, 



for all i Pl ,ip 2 ,i P3 GV([0,T) x Si). 

Remark 1. Note that we can easily check that Definition 1 1 . 1 1 makes sense. Further- 
more, observe that Definition 11.11 implies that dtUi belongs to L 2 (0,T; (H 1 (S})Y), 
so that Ui G C{[0,T};L 2 (Slj) for i = 1,2,3. 

Remark 2. A classical way to prove the existence of weak solutions in the sense of 
(jl.6|l - (|1.8j) . is to use Faedo-Galerkin method like the system studied in [I] or in [14] . 
On the other hand, the proof here, of convergence of the numerical scheme, implies 
the existence of weak solutions of (|l.ljl - (|1.3|) . Additionally, a proof of uniqueness 
of the weak solution is given in the appendix. 

Following [9] , we now give a precise definition of the finite volume scheme for the 
nonlocal SIR model. Let SI be an open bounded polygonal connected subset of M 3 
with boundary dSl. Let Sin be an admissible mesh of the domain SI consisting of 
open and convex polygons called control volumes with maximum size (diameter) h. 
For all K G SIr, let by xk denote the center of K , N(K) the set of the neighbors 
of K i.e. the set of cells of SIr which have a common interface with K, by Ni nt (K) 
the set of the neighbors of K located in the interior of SIr, by N cxt (K) the set of 
edges of K on the boundary dSl. Furthermore, for all L G N(K) denote by d{K, L) 
the distance between Xk and xl, by (Jk,l the interface between K and L, by rj^ £ 
the unit normal vector to <jk,L outward to K. For all K G SIr, we denote by m{K) 
the measure of K. The admissibility of SIr implies that SI = Uk£Q r K, K n L = 
if K, L G SIr, and K ^ L, and there exist a finite sequence of points (xk)k£Q r 
and the straight line xkxl is orthogonal to the edge ctk,L- We also need some 
regularity on the mesh: 

d(K,L) 
mm — - > a 

K£Ci R ,LeN(K) diam(A) 

for some a > 0. 

We denote by Hh(Sl) C L 2 (Sl) the space of functions which are piecewise constant 
on each control volume K G SIr. For all Uh G Hh{Sl) and for all K G SIr, we 
denote by uk the constant value of Uh in K. For {uh,Vh) G {Hh(Si)) 2 , we define 
the following inner product: 

/ \ 1 m ( a K,L) , w v 

Kefin LeN(K) 1 ' ; 

corresponding to Neumann boundary conditions. We define a norm in Hh(Sl) by 

\\ u h\\ Hh (Q) = (( u h,U h ) Hh ) 1/2 . 

Finally, we define Lh{Sl) C L 2 (Sl) the space of functions which are piecewise con- 
stant on each control volume K G SIr with the associated norm 

{Uh,V h ) Lh{n) = Yl m ( K ) u KV K , W u h\\l h (Q) = m ( K )\ u K\ 2 , 

Ken R Ken R 
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for (uh,v h ) £ {L h (Q)) 2 . 

Next, we let K E Qr and L E N(K) with common vertexes (ae,K,L)i<£<i with 
I EN*. Next let T k ,l (respectively T™^ for a E N cxt (K)) be the open and convex 
polygon with vertexes (xk,xl) (xk respectively) and (at t K,L)i<t<i- Observe that 



- UKen R I {^LeNiK^K.Lj^^JaeN^tiK^K^ 



The approximation V^u/i of Vit is defined by 




I m((TK J.) r 



o ifaser^. 



for all if € f2fl. 

The next goal is to discretize the problem (|1.1[) - (|1.3[) . We denote by V an 
admissible discretization of Qt, which consists of an admissible mesh of f2, a time 
step At > 0, and a positive number N chosen as the smallest integer such that 
NAt > T. We set t™ = nAt for n E {0, ... , N}. 

We approximate the nonlocal SIR model in the following way: Determine vectors 
(u™ K )Ken R for i=l,2, 3, such that for all K E £Ir and n E {0, . . . , N — 1} 

1 



(1.9) U ^ = m(K) I u i;o{x)dx, i= 1,2,3, 



A' 



(1.10) A GO ft LeN(K) y ' ' 

n+l n / \ 

(1.11) At A oG Q h LEN(K) a(K ' L) 



m(K) (v(uZ+,u n 2 + 1 > + ,u n 3 + 1 ' + ) ( 7 + M K+ X ) = 0, 



f T S\ U 3,K / n \ m \ (J KX) ( n+1 n+l\ 

M121 W( } At H2>3,a J £ -d(KjS) {u ^ " U 3,a) 

- m(if ) 7 uS )A . = 0. 

To simplify the notation, it will always be understood that when h is sent to zero 
then so is At, thereby assuming (without loss of generality) a functional relationship 
between the spatial and temporal discretization parameters. This is not a real 
restriction on the time step, but it allows us to write instead of At" f° r 

i = 1, 2, 3 , "h — > 0" instead of "/i, At — > 0", and so forth. 

For the sake of analysis, we introduce the "piecewise constant" functions 

(1.13) u hh (t,x) =u?+\ i = 1,2,3, 

for all (t, x) E (nAt, (n + l)At] x K, with if e SIr and n 6 {0, . . . , N - 1}. To 
simplify the notation, let us write for the vector (wi^, v,2,h, U3,h)- Our main 
result is 



FINITE VOLUME SCHEMES FOR NONLOCAL REACTION-DIFFUSION SYSTEMS 5 



Theorem 1.1. Assume Uifl G L 2 (fl) for z = 1,2,3. Then the finite volume solution 
Uh, generated by (|1.9I) and (|1.10|) - (|1.12l) . converges along a subsequence to u = 
(1*1,1*2,^3) as h — > 0, where u is a weak solution of (|1.1[) - (|1.3[) . The convergence 
is understood in the following sense: 

Ui : h Ui strongly in L 2 (Qt) and a.e. in Qt, 

^hUi,h -> Vu, weakly in (L 2 (Q T )) 3 , 

v(ui,h, U 2 ,h, u 3 ,h) -> ct(mi, it 2 , u 3 ) strongly in L 1 (Q T ), 

for 1 = 1,2,3. 

The remaining part of this paper is organized as follows. The proof of Theo- 
rem [O] is divided into Section [2] (existence of the scheme), Section |3] (basic a priori 
estimates), Section |4] (space and time translation estimates), and Section [5] (conver- 
gence to a weak solution). In section [6] we give some numerical examples. Finally 
in Appendix we prove the uniqueness of the solution using duality techniques. 

2. Existence of the finite volume scheme 

The existence of a solution to the finite volume scheme ()1.9|) - (|1.12p will be ob- 
tained with the help of the following lemma proved in [12j and |16j. 

Lemma 2.1. Let A be a finite dimensional Hilbert space with scalar product [■, ■] 
and norm \\-\\, and let V be a continuous mapping from A into itself such that 

[P(0,$>0for U\\ =r>0. 

Then there exists £ 6 A with ||£|| < r such that 

v(0 = 0. 

The existence for the finite volume scheme is given in 

Proposition 2.2. Let T> be an admissible discretization of Qt- Then the problem 
(jl.9p - (|1.12p admits at least one solution {u™ K ,uV^ K ,u^ x)(K,n)en a x{o,....N}- 

Proof. First we introduce the Hilbert spaces 

E h = {H h (n)nL h (n)f, 

under the norm 

3 3 

\\ u h\\% h : =Hl! u ^Hk(fi) + E E m(K)\u hK \ 2 , 
i=i i=i if£0 H 

where u h = u 2 ,h, u 3 ,h)- Let $ h = (tpi,h, f2,h, V>3,h) G E h and define the 

discrete bilinear forms 

3 

i=l 

&, 1 K+ 1 ,$, l )= E m(ir)(a(<+ 1 '^<+ 1 '^«^ + )+ M <+ 1 ) ¥)1 , x 

EIts\( I + n+l,+ n+l.+ \ / . \ n+l\ 
m(K) {o-iu^x, u 2 K , u 3 K ) - (7 + n)u 2 ^ K J ip 2 ,K, 

Ken R 
Ken R 
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and 



E <Ko 



a h (u n h +L ,<f h 



2 

i=l 



E E ! ?7^r«,i 1 -<i 1 )«i 1 -<i^ 



Kens lgn(k) 



d(K,L) 



Multiplying (jl.lpp . (jl.lip and (|1 . 12[) by ^.r-, <P2,k, <p3,K> respectively, we get the 
equation 



= 0. 



Now we apply the Lemma [2TT1 for proving the existence of u^ +1 for all K 6 and 
n£ {0, . . . , iV}. We define the mapping V from into itself 



+ a h (u n h + \$ h ) + b h (u n h +1 ,$ h ), 

for all S Eh. Note that it is easy to obtain from the discrete Holder inequality 
the following bounds: 

a h (u h ,v h ) < C\\u h \\ Eh \\v h \\ Eh , 
T h (u h ,v h )<C\\u h \\ Eh \\v h \\ Eh , 
b h (u h ,v h ) < C\\u h \\ Eh \\v h \\ Eh , 

for all Uh and in Eh- This implies that ah, Th and bh are continuous. The 
continuity of the mapping V follows from the continuity of the discrete forms dh{-, ■)> 
T h (; •) and b h (-, •)■ 

Our goal now is to show that 



(2.1) 



[p( U r i ),uri>o far k 



n+1 1 



r > 0, 



for a sufficiently large r. We observe that 



(2.2) 



b h(U h ,U h ) 



,n+l 
L i,K 



+ a, 1 K +1 ) < +1 ) 



^E E w^ 1 - 

i=i k en R 
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It follows that from the definition of a and (|2.2| and Young's inequality that 



3 3 

^E E m(K)\ulP\ 2 + J2 M " 



=1 Ken R 



,71+1 



2 

i? h (n) 



n+l,+ n+l, + \ 7i+l, + , 

, u 2if , u 3if )u 1K + /i 



Ken 



,71 + 1 
^l.if 



m{K)(a{u n ^ 

:en R 

Ens\( I n,+ 77+1. + 71+1, + \ , \ 7i+l\ n+1 

m{K)\a{u x ' K ,u 2K ,u 3J< ' ) - (7 + fi)u 2 ^ K ju 2 ^ K , 

" fin 



3 



i=i Ken n 



7=1 KGfij- 



-4- V m(iH 
8Ai <^ 1 ; 



A7 E 



8Ai 



71+1 



77 + 1 

i 2,X 



77 + 1 



i2 

l-fffc(fi) 



7=1 ffGfiij 



77,+ 



,77+1 

X 7,A- 



7=1 

-d(At,a) ^ m(if) 
Jfe!i H 

-C 2 (Af, 7 ) £ ™(^)K*| 2 
3 



7=1 jYgSIj 



This implies that 
(2.3) 

3 



>min{ — ,M^M 2 ,M z }\\nl +1 f Eh 

3 

-2max{c 1 (At,a),C 2 (At ) 7) ) C 3 (At)}^ ]T m(if) 



7=1 ife^j 



Finally, for given u" h , u 2h and h , we deduce from (|2.3j) that (|2.1|) holds for r 

= r) . Hence, we obtain the existence of at least 



77+1 1 



large enough (recall that ||ujj 
one solution to the scheme (|1.9[) - (|1.12p . 



□ 



2.1. Nonnegativity. We have the following lemma. 

Lemma 2.3. Let (u™ K , u 2 K , U3 K)Ken R ,ne{o,...,N} be a solution of the finite vol- 
ume scheme ljl.9j) . (jl.lOp . (jl.lip and (|1 .12j> . Then, (u™ K , u 2K , u 3 K )K^n R ,n^{o,...,N} 
is nonnegative. 
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Proof. Multiplying (flTO]) by -Atu^' , we find that 
(2.4) 



m(K)u^-(u^-ul K ) + a 1 ( JT At £ ^A {u ^ - u^)u^ 



K en h 



LeN(K) 



fTS\\j-( I n+l.+ 7i+l,+ n+l,+ \ n+l _ , ri+l\ 

m(i^)A^cr(u 1K ,u 2 y ,u 3J( ' +M«i,kJ u i,k 



We know that u^-' ' = u 
With this, we deduce 

(2.5) 



n+l _ „,n+l+ „,n+l 



A" 



— and (a+ - 6+) (a" - &") < for all a, & £ 



JV-l 



n=0 Aefij? LeJV(A) v ' y 
( £ U ™ A ) w _! 



m{o- K ,L) , w+ i _ „ + i^ n+l- 
V"l,L U 1,A/ U 1,A 



K en h 



ai | 



A en h 



> 



"'i E M "^o 

A" en h 



E Ai E E Ki x -<£)Kr 



1,A J 



„=0 Jf£S! H LeN(K) 



JV-l 



E A 'E E K 

n=0 AeO R LeJV(A) 



n+l~ _ „,n+l~ 
L "l,A 



> 0, 



and 



(2.6) 



JV-l 



E A * E 



n+l,+ ,.n+l,+ ..n+l,+\ _j_ ,,„,n+l \„,n+l~ 



n=0 Aefin 



A ' "2, A ' "3, A 
JV-l 



- M e A * E 

n=0 AeSl R 



n+1" 
*1,K 



< 0. 



Let / S C 2 function. By using a Taylor expansion we find 



(2.7) 



/(&) = /(a) + /»(& - a) + -/"(£)(& - «) 2 , 



for some £ between a and 6. Using the Taylor expansion (|2.7p on the sequence 

f(u^ + jl) with /(p) = f sds,a = and & = v% K . We find 
Jo 



- |2 



(2.8) « 



7t+l /..H+l 



n.A 



,n+l- 



A "%AJ 



2 2 2 

We observe from the definition of / that f"(p) > 0, which implies 

2 

(2.9) <+ x «£-<k)< 



A 



n,A 
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Now, using (|275 [) -pT9 j) to deduce from §2A$ 



N-l U, n +S 



"1,K 



E 

(2.10) n=0 



This implies that 
(2.11) 



"1,K 



2 

N-l 



n=0 



+ M E A * £ m(JQ 



n=0 KESlf 



,n+r 

4 1,K 



"1,K 



< 0. 



< 0. 



Note that (|2. 12[) is also true if we replace N by no € {1, ...,N}, so we have 
established 



(2.12) 



l l.K 



< 



l l,K 



Since u\ K is nonnegative, the result is vlf^ for all < n < N — 1 and all K G 
On the other hand, multiplying (|l.llj) by —Atu^^ , and along the same lines as 



z i if ' we obtain the nonnegativity of discrete solutions for all < n < N — 1 
and all K £ Finally, the nonnegativity of w^^ 1 , is given by a maximum princi- 
ple proved in [5]. In fact the reaction terms m(K) 7 K of the equation (|1.12p do 
not depend on and it is a nonnegative term. We assume the nonnegativity of 
( u i Ki u 2 K' u 3 k) ano - we a PPly the discrete maximum principle for the third equa- 
tion (|1.12[) in order to prove the nonnegativity of . Then, using an induction 
on n yields, we conclude. □ 



3. A PRIORI ESTIMATES 

The goal is to establish several a priori (discrete energy) estimates for the finite 
volume scheme, which eventually will imply the desired convergence results. 

Proposition 3.1. Let (uf K )^ e Q R . n e{o,....N}> i — 1)2,3, be a solution of the finite 
volume scheme (|1.9p - (ll.l2p . Then there exist constants Ci, 62,63 > 0, depending 
on Q, T, Ui_o and a such that 



(3.1) 



max m{K) \u™ K \ 2 < C\, 

l£{l,...,N} ' 



Ken T , 



(3.2) 
and 



N-l 



J E At E E 



71=0 



K£Sl R L^N(K) 



d(K, L) 



n+1 

L i,K 



,n+l 

l i,L 



<c 2 , 



N-l 



(3.3) J2 At E m{K)[a{ul K ,u^,u^) + a(<+ 1 ,^+ 1 ,<+ 1 ; 
fori =1,2,3. 



<C,. 
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Proof. We multiply (fTTOll . (fTTTTjl and (fl~T2l) by Atu™+ 1 , Aiu^ 1 and Atw^+ 1 , 
respectively, and add together the outcomes. Summing the resulting equation over 
K and n yields 



E\ + E2 + £3 — -E4 



where 



3 JV-l 
i=l n=0 if G fin 



^ m(K)(<+ 1 -< K )<+ 1 , 



3 JV-l 



i = -EEM E E 



i=l n =0 



K ^n h Ken R lgn{k) v ' ' 



m(a K ,L) , n +l _ n+lN n+l 
\ a i,L a i,K ) a i,K > 



£3 



AT-1 

E 



K£fl R \ 



+ (7 + M) 



,71+1 

x 2,if 



iV-1 



^ = E At E "•(^(^K^^-^K^^ 

n=0 Jfefia 

From the inequality "a(a — b) > \{a 2 — b 2 )" , we obtain 



„7i+l 

if u 3.if 



3 N-l 



£ i = EE E ^(^-W 

i—l n =0 if Sfi H 
-. 3 N-l 

>Jee E 

1=1 n=0 ifeSln 

= ^E E 



n+l 



l N I 

i=i ifen B 
Gathering by edges, we obtain 
3 w-i 

^ = -EEM E <k ) e e 

j=l n=0 K a en h Ken R LeN(K) 

m(<TK,L) 



H,K\ 



m{o K ,L) , n+ l _ n+ u n+ i 

d(K,L) [U *> L ' K 



3 M- JV " 1 

^EtE a * E E 



z=l n=0 Ken R L<EN(K) y ' ; 



Observe that from nonnegativity of {u^ k )k en R ,ne{o,...,N} f° r * = 1> 2, 3, we get 

£3 > 0. 
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We use Young's inequality to deduce 



JV-l 

n =o Ken™ 

N-l 



=0 A'esij 



JV-l 



l 2,K 
N-l 



N-l 



<aJ2^J2 m(K)\u^\ + ~£At £ 



n=0 AgQr 



|E A * E ™w 



n=o agOr 



..n+l 
1 3,K 



<( a + |)]T At £ m(K) 



=0 AGO* 



..n+l 
1 2,K 



7 



At J2 m ( K )\ u li<\ 



Ken, 



N-l 



+ 1 E A< E m ( K ) 

n =a Ken R 
Collecting the previous inequalities we obtain 

^E E m^d^i 2 -!^! 2 ) 

i=i Ken R 

3 , r N-l 



..n+l 
1 3,K 



+ 



EfE^E E 



(3.4) 



m(a KtL ) 

— 2 ^ ~" ^ ^ d(K, L) 
i=i n =o Ken R LeN(K) v ' ' 



a 



n+l _ ..n+V 
,A a i,L . 



N-l 



<(a + |) E Ai E + \ At E m ^ K* 



n=0 AgOr 



N-l 



J E A * £ 



n=0 Aefifl 



Aerij 



n+l 
^3 : A 



which implies 



£ E -wKa| 



i=l AgSIr 



3 JV-l 

<£ £ m(if)|<^| 2 + (a+ At ^ 

i=l AgOr n=0 AgSIr 

JV-l 

+ |At J2 ™(*0Ka| 2 +|£ a * E ™w 



(3.5) 



Ken R 



71=0 Aesi f 



,n+l 
*2,A 



,.T»+1 

X 3,A 



Clearly, 



E m ( K )\ u ii<\ 2 <\\^,o\\ 2 L 2(n) fori = 1,2,3. 



AgOr 



In view 



-LV R. 

iew of ()3.5|) . this implies that there exist constants C4, C5 > such that 



(3.6) ]T ^ m(JT) 

i=l AgOr 



<k\ 



3 JV-l 

< C 4 + C 5 ^ J] At m ( K ) 
i=i n =o Ken R 



..n+l 
I i.K 
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Note that (|3.6|) is also true if we replace N by no G {1, . . . , N}, so we have 
established 

3 2 3 no — 1 n 

(3.7) £ ]T m(JO u& <C7 4 + C 5 £ £ Ai £ m(X) u "+? 

i=l KeH fl i=l n=0 Keil a 

By the discrete Gronwall inequality (see e.g. [ID]), we obtain from (|3.7j) 

3 



(3.8) 



E E »»(*) 



=i Keu h 



l i.K 



<C 6 , 



for any no G {1, . . . , AT} and some constant Ce > 0. Then 



max m(K)\uf K \ < C 6 



i=i Ken 

Moreover, we obtain from (|3.4[) and (|3.8p the existence of a constant CV > such 
that 

3 JV-1 



EE A * E E 



m((JK,L) 
d(K, L) 



i=l n=0 Ken R L£N(K) 

Finally a consequence of (| 1 .4|) and (|3.ip is that 

||o-(wi,h)«3,/i.«3,fe)lli2(Q T ) < C 8 , 

for some constant C% > 0. 



< C 7 . 



□ 



4. Space and time translation estimates 

In this section we derive estimates on differences of space and time translates of 
the function Vh which imply that the sequence Vh is relatively compact in L 2 (Qt)- 

Lemma 4.1. There exists a constant C > depending on fl, T , 1*1,0, M 2.o, ^3,0, 
a, 7 and /1 such that 



(4.1) 



K, h (t, x + y) - tii.fcft x)\ 2 dxdt<C \y\ (\y\ + 2h), i = 1, 2, 3, 



/n'x(o,T) 

/or aiZ 1/el 3 iwzf/i !!' = {1 6 0, [a:, a; + y] C 0}, and 



(4.2) 

'Slx(OJ-r) 

for all t G (0,T). 



+ t, x) — Ui : h(t, x)\ 2 dx dt < C(t + At), i = 1, 2, 3, 



Proof. The proof is similar to that found in, e.g, [5]. 

Proo/ 0/ ijlJj) . Let y G K 3 , x G fi', and L G AT(iT). We set 



Next, the value 



(4.3) 



1, if the line segment [x,x + y] intersects ctk,l, K and L, 
0, otherwise. 

U 

c aK L is defined by c„ K L — — - ■ r\x l with c aK , > 0. Observe that 

\y\ 



Xvk,l( x ) dx < m(a K ,L) I y\ c <? 
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Using this, we obtain 

\ui, h (t,x + y) -ui lh (t,x)\ < ^2 X<t K ,l{ x ) 



n+l _ n+l 
U 1.L U 1,K 



To simplify the notation, we write 

instead of 



E 



0jc,l {(K,L)£n 2 R ,K=£L,m(a K:L )^0} 

By the Cauchy-Schwarz inequality, we get 

\ui, h (t,x + y) -uiM{t,x)\ 2 < ^2 Xa K , L (x)c aK L d(K,L) 

<?K,L 

(4.4) i - ■ - 2 



E 



\v n+1 - v n+1 
"l,L "l.A 



~X&K,L ( X )- 



Note that 
(4.5) 



E X° K , L ( x ) c °K, L d(K,L) < \y\ + 2h. 



Using (|4~3ll . (|P|) . and (|4~5ll . we deduce 



|tti,/t(t,a; + y) - 



(4.6) 



<(|y|+2^)^ A*£ 

n=0 ok,l 
N~l 

<\y\(\y\ + 2h)J2^J2 



N-l 



n+l _ n+l 
"l,L "l.if 



n=0 ctk.l 



c aKL d{K,L 

m((j K . L ) 
d(K, L) 



Xcr KiL (x)dx 



,,«+! _ ,,"+1 

U 1,L U 1,K 



Then, from (pH)) and (HI))) , we deduce P~Tj) . 

PVoo/ o/ g3). Let r G (0, T) and t e (0,T- r). We have 



B(t) = \ui t h{t + r, x) -ui t h(t,x)\ dx. 
Jn 

Set n (t) = E(t/At) and m(f) = £((t + r)/Af), where 



= 7i for x e [n, n + 1), n E N. 



We get 



B(t)= ^ m(JQ 



U 1,A "l,K 



This implies 



(t) _ n (t) 
A' U 1,A 



Lefl a 



53 m(^)K+ 1 -< )X ). 

i<(n+l)Ai<f+T 
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Using the scheme (|1.10p . wc obtain 



B{t) = Ai E K 1 * ^ <t 

t<(n+l)At<t+T Ken R 



n+1 _ n+1 
"l.L U 1,K 



(4.7) x h( E «?.*.) E d(x £) 

- mm (<7K+ 1 ,<+ 1 ,^+ 1 ) + 

We observe that we can rewrite (14.71) as 

B(t) =i E ME »?,«„)£ E 

t<(n+l)At<t+r K ett h K£fl R LeN(K) v ' ' 

+ E A * E 0^ - 

i<(n+l)Ai<t+T Ken n 

x |-m(^) <+\ + j . 

We use the basic inequality "ab < ^a 2 + ^b 2 " to deduce 

B(t) < i(B!(t) + ^B 2 (t) + iB a (t))+B 4 (*), 

with 

Bl «)= e E 

i<(n+l)At<i+r A"o60 h Ken R L£N(K) v ' ' 



*1,K "■!,£ 



t<(n+l)At<t+r ifo6^h Ken R L£N{K) y ' 



m(<7 K ,L) 

d(K, L) 
m(<7 K ,L) 



«.(*)= E ME < Ko )E £ 

i<(n+l)At<i+r if en h KeQ R LeN{K) v ' ' 

B 4 (t) = - £ At ™(K)Wl 

t<(n+l)At<t+r Ken R 



2 

,«l(t) _ nl(t) 

X 1,K a l,L 



2 

no(i) _ »o(t) 
l l,K a l,L 



m(*) _ „no(t)\ 



x(aK+\<+\<+ 1 )+^<+- 1 ). 

Next, we introduce the characteristic function \ defined by x(n, ii, t 2 ) = 1 if 
t\ < (n + l)Ai < ^2 and x(^>ti>t2) = otherwise. Then we have for any sequence 
(a")„ e N of nonnegative numbers that 

(4.8) / J2 a n dt< J2 °" / x(n,t,t + r)dt <t a " 

t<(n+l)Ai<t+r n=0 J n=0 
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and for any £ G [0, r] 

r T~r 



(4.9) 



a £((t+o/At) 

t<(n+l)At<t+T 
E (-Kt) p(rn+l)At 

^ E / a? ™ E xCM-fct-f + r)* 

m=0 JmAt „ =0 
B (st) »Ai B (st) 

= E / a ™ E x(»-"».*-£.*-£ + T)<tt 

m=0 ^° n=0 

= E E / a m x(-™,t-e,t-e + r)dt 

m= n=0 "'-"A* 

= E a ™ / x(-m,*-e.t-f + r)dt 

m=0 ^ R 

B (3t) ,.f-(n+l)At B (st) 

= am dt = T E a "- 

m=0 ^+r-(n+X)At n=0 



From (14. 8|) . we deduce 
-t-t JV-i /.T-T 



[ B 1 (t) dt<J2 At [ X(n,t,t + T ) ai (j2 u Ik 

E E 



n=0 u ifoGO,, 



K&Q R LGN(K) v y 



U 1,K "l,L 



n=0 if e£l h (K,L)£n 2 H 



d(K, L) 



n+l _ n+l 



Using ()3.2|) . this implies that there exists a constant Cio > such that 

r T-T 

"^10- 



,T-r 

/ Bl(t)dt<TCl, 
./0 



Now we consider B 2 (t) and B 3 (t). We use (gUJ) with £ = r for S 2 (i) and g^J) for 
^(t) to obtain 

/ B 2 (t) dt < r At a x ( ]T «J jJfo ) £ £ 



m(a K ,L) n+l ..n+l 



u l,_fS" U 1,L 



and 

„T_ r A - I 



m(ffA-,L) 

d(K,L) 

7i=0 -fS"o£f2h K£U R LeN(K) v 7 



/ %(*) A < r X) At O! ( 5] «U„) E E 



"l,K "l,L 



In view of (|3.2[) , we deduce 

f>T-r f T- 



/ B 2 (t)dt<TC 11 , f B 3 (t)dt<rC 12 , 
Jo Jo 



for some constants Cn, C12 > 0. 
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Next, we intergrate Z? 4 (i) from to T — r. Using the inequalities "ab < ^ + 



and "(a - b) 2 < 2a 2 + 2b 2 '", we obtain 



n+ 
K 



t<(n+l)At<t+7 



Ken R 



m(t) 



^ E A< E m ( K h 'm.a 

t<(n+i)At<t+T ifenfl 
From and ([XT]) , we have 

£ Ai £ m(K) 

t<(n+l)At<t+r KeVln 

for some constant Ci 3 > 0. We use (|3.ip to deduce 



710 (*) 



+ C(a,/i) 



,71+1- 
4 1,K - 



T» (t) 



eft < tC 



13) 



£ Ai £ m(tf) 

t<(n+l)Ai<t+r Jfefis 

for some constant Cm > 0. Hence 



.71+1 



< tC u , 



B 4 (t) <tCi 5 , 

for some constant C15 > 0. Reasoning along the same lines for u\ t h yield (14. ip and 
(|4.2p for U2 h and U3 ^- This concludes the proof of the lemma. □ 



5. Convergence of the finite volume scheme 



The next lemma is a consequence of Lemma 14. H and Kolmogorov's compactness 
criterion (see, e.g., [2J, Theorem IV. 25). 

Lemma 5.1. There exists a subsequence of = (iti,/,, M2./1, M3./1), ^of relabeled, 
such that, as h — + 0, 

(V — * strongly in L 2 (Qt) and a.e. in Qx, 
(5.1) fiij VhUi,h — * Vwj weakly in (L 2 (Q T )) 3 , 

(Hi) o-(ui,h,U2,h,U3,h) — > cr(ui,u 2 ,U3) strongly in L 1 {Qt), 
fori =1,2,3. 

Proof. The claim (i) in (|5.ip follows from Lemma 14.11 and Kolmogorov's compact- 
ness criterion (see, e.g., [2], Theorem IV. 25). The proof of the claim (ii) will be 
omitted since it is similar to that of Lemma 4.4 in [3], we refer to the proof of this 
lemma for more details. The claim (iii) follows from Vitali theorem. □ 

Our final goal is to prove that the limit functions ui,U2, U3 constructed in Lemma 
15.11 constitute a weak solution of the nonlocal system (|l.ip - (|1.3p . 

We start by verifying (| 1 . 6[) . Let T be a fixed positive constant and ipi G 
X>([0,T) x fi). We multiply the discrete equation ([LTD]) by At(p 1 (t n , x K ) for all 
K E and n € {0, . . . , N}. Summing the result over K and n yields 

T 1 +T 2 + T 3 = 0, 
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where 

N-l 

n =o Ken R 

N-l , . 

t 2 =-y: A tai ( £ ui Ko ) y: e Sttki 1 - 

n=0 K en h K£fl R LeN(K) 1 ' ' 

T 3 = E Ai E ™(^)( CT K+\<+\<+ i )+ M <+ i )^ 1 (r ;a;i ,). 

n=0 KGO K 

Doing integration- by-parts, keeping in mind that <pi(T, xk) = for all K G f2#, 
we obtain 

JV-l 

Ti = -E E "i^Xi 1 ^!^ 1 ,^)-^!^",^))- E rn{K)ul K yi{Q,x K ) 
n =o Ken R i<en R 

N-l ~t n+1 p r 

= ~E E / / u ijidt<Pi(t,x K )dxdt - ^2 / "i,o(£Vi(0>a ; K)<fo 

='■ —T\ t i — T\_2- 
Let us also introduce 

^ = -E E / u r {+ K 1 d t vi(t,x)dxdt- u 1;0 (x)ip 1 (0,x)dx 

n=o Ken R Jk •' n 

Then 

T , i,2-T' 1 *2= E / ui^(x){tfi(Q,x K ) - yi(0, a;)) da;. 

Mr''' 1 ' 

From the regularity of there exists a positive constant C such that 

|(pi(0, a:^) - Vi(0, £)| < C ft,, 

which implies 



\Ti.2 - T* 2 \ < Ch 2~] / \uifl(x)\ dx. 



(5.2) 

Sending h —> in (|5.2p . we get 



lim Ti 2 — 7T 9 = 0. 
?i->q 1 ' ' 1 

Observe that 
^1,1 - 



E E I / / u^dtVi{t,x K )dxdt- f f u^d t <pi(t,x)dxdt 

N-l t n+l 

E E u Tk J J [d m {t,XK)-d m (t,x))dxdt. 
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Using the regularity of dt^pi and Holder's inequality, we obtain 

/N-i „ \ 1 / 2 

\T 1A -T* lA \ <C(h) At E m W 
\n=o Ken R 

where C(h) > is a function satisfying C(h) — > as h — > 0. From (|3.ip we deduce 











U 1,#T 





lim Ti i -T* = 0. 



Next, we define Id and T 2 * by 
Id 



/ a i[ u\dx) \7ui ■ Vifi dx dt, 
Jo vn ' in 

^2* = / a i( / "i./i^I / Ui^Atpi dx dt, 
Jo vsj ' Jn 



/ a i ( / uidx) uiAipidxdt. 
Jo vn ' Jq 



Integration-by-parts yields 
Id 

On the other hand, using the convergence results of Lemma l5.il and taking into ac- 
count the assumption (jl.5|) . it is easy to prove that a\ (J„ u i.h dx) — > ai (Jq mi da:) 
strongly in £ 2 (0, T) and J n Mi^.A^i da; — » — J n Vui • V^i da; weakly in £ 2 (0, T), as 
h — ► 0. Thus, there holds 

T 2 * -> -7 D as h ->■ 0. 

Note that 

,Ld7 

/(if) 



= E fl i( E E E / / v ^ ■ 

n=0 X 6n h K&n R LGN(K) " *" " &K,L 

= -2E fll (E "i,KdJ E E ( u i% 1 - u ijc) / V<Pi-VK,Ld-y 



n=0 KoGn,, K£Sl R LeN(K) 



and 



iV-l „ n + 1 - 1 ,™+ 1 



T 2 = -^A tal ( £ < Ko ) £ E X t K l< K <Pi(t n ,x K ) 

n =o K en h Ken R LeN(K) ^ ' ' 

N-l , , 



n=o /f eo ft Ken R LeN(K) v ' ; 



x «r - <^)(^(r,a; L ) - <pi(t n , xk)). 



Hence 



r 2 +^=^ fl i( E u U„) E E m^^Ki 1 -^ 1 ) 

n=o i<" ef2fe Ken R LeN(K) 

» cpi(t n ,*L)-<Pl(t n ,XK) 

d(K, L) 

V<pi ■ r\ K ,hd^ . 



m{a KyL ) J t 



&K.L 
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Since the straight line {xk,xl) is orthogonal to ok.l, we have 

xk - XL = d(K, V)t)k,l- 
This implies from the regularity of ipi that 

(pi(t n ,x L ) - (px(t n ,x K ) 



and so 



(5.3) 



d(K, L) 

Vl {t n ,x L )~ Vl {t n ,x K ) 
d(K, L) 



Vipi(t n , x) ■ rjK,L with x between xk and xl, 



1 



m(a K ,L) J t 



< CAt h, 



for some constant C > 0. Using (|5.3p and (|3.3p . we deduce 



lim T2 = — / a± ( / i*i da;^ / ui A</?i dx dt = / a i ( / u i / ^Ju\-Wtp\dxdt. 

fi(ui, U2, U3) + [lUijtpi dxdt. 



Now, we show that 



lim T3 

For this purpose, we introduce 
N—i 



eft 



and 



n=0 Kefi R 



N-l 



T * -EE 

n=0 KgQr 



N-l 



/i(u"^ — ui)(pi(t, x) dxdt. 



dx dt 



+ E E 

We have for all x G X and £ S [i n , i n+1 ] that 

(5.4) \<Pi(t n ,x K )-<px{t,x)\ <C(At + h), 

and thus, thanks to (|3.ip and (|3.3p . 

N-l 

|T 3 ,i| < C(At + h) J2 At m(K)(a(u^ K \u^\u^ K 1 ) + ^+ K })<C(At + h). 

n=0 Ke(l a 

Hence, T3.1 — > as h — > 0. We also have 



\Ts,2\ < C / / (\a(ui t h,U2 t h,U3,h) - o-(ui,U2,m)\ + \ui t h - Ul\ ) 



dx dt. 



10 Jci 

Therefore from this and (15.11). we deduce 



|T3 5 2 1 tends to zero as h — > 0. 

This concludes the proof of (|1.6p . Reasoning along the same lines as above, we 
conclude that also (fTTTl) and (Oil hold. 
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6. Numerical Examples 

6.1. Example 1. SIR model simulations. In this section we consider a sample 
square domain = (0, 1) x (0, 1) and we show the behavior of the solution for 
different models of nonlocal functions <Zj, i = 1, 2, 3. 




FIGURE 1. SIR model with oi = a 2 = a 3 = 1/10: (a) Beginning 
infected population (t = 0.025 sec); (b) Susceptible population 
(t = 0.5 sec); (c) Infected population (t — 0.5 sec); (d) Recovery 
population (t = 0.5 sec). 



We consider here a uniform mesh given by a Cartesian grid with N x x N y control 
volumes and choosing N x = N y = 300 for all simulations. Obviously, it is possible 
to consider also unstructured meshes, but we will confine here to uniform mesh for 
simplicity of the simulated models. The discretization in time is given by Nt = 100 
time steps for T = 0.5. That is, St = T/N t and m(K) = l/(N x N y ). The parameter 
of the SIR model are given by a = 2.0, fi = 0.01 and 7 = 1.0. 
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Figure 2. SIR model with <n = a 2 = a 3 = s/10: (a) Begin- 
ning infected population; (b) Susceptible population; (c) Infected 
population; (d) Recovery population. 

The initial condition are given by 

{ui,o(x,y) = e ; 
U2,q( x iV) = B^sech(/3(x - Xj))sech(/3(y - %)); 
3=1 
u3,o(x,y) = 0.0. 

withe = 0.01, B = 5000, j3 = 2000, = (0.25.0.25), (x 2 ,y 2 ) = (0.125.0.125), 

(x 3 ,y 3 ) = (0.125.0.375), (x 4 ,y 4 ) = (0.375.0.125) and (x 5 ,y 5 ) = (0.375.0.375). 
These initial conditions represent one hand to a susceptible population initially with 
a low density and constant throughout the domain. On the other hand, 5 pockets of 
high density infected population are located in the quadrant [0, 1/2] x [0, 1/2] which 
will diffuse the epidemic desease on the rest of the domain. Finally, we assume that 
there is no initially presence of recovery population. 

In a first simulation we consider a model which the diffusion rates do not depend 
on the population, that is a^, i = 1,2, 3 are constant and equal to 0.1 (see Figure 
[1}. On the other hand, in figure [2] we observe another simulations with a nonlocal 
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Figure 3. Example 1: evolution in time of the diffusion terms 
a, (J Ui dx^j , with i — 1, 2, 3. 



diffusion where the three diffusion rates are linearly proportional to the total mass 
of each population (&i(s) = 02 (s) = 0,3(3) — 0.1 s). Outside the nonlocal diffusion, 
we consider the same parameters for both simulations. We remark that second 
simulation with nonlocal diffusion violate the assumption (|1.5p . In fact, they cor- 
respond to degenerated parabolic cases. In order to ensure the convergence of our 
numerical example, we replace in a practical way, the diffusion rate coefficients by 



(6.1) 



M if s > M 
Oj(s) = { s if e < s < M 

£ if S < E, 



with M = 10 4 and e = 10~ 4 . Figures [T] and [5] represent the simulation at time 
t = 0.025 for the localized infected population (see picture (a) of both figures) and 
the simulation at time t = 0.5 for the different populations (see pictures (b), (c) 
and (d)). We remark that different behaviours occur for the susceptible population, 
observing a first extinsion of this population in a different zones of the square 
[0, 1] x [0, 1] for each case. Finally, we observe in Figure El the evolution in time 

of the diffusion (^J U{ dx^j , with i = 1, 2, 3 for the simulation with the nonlocal 

diffusion. 

Unfortunately, it is not possible to observe patterns formation for this SIR- 
diffusion model like some results obtained in [T3]. That is because the infected 
population have his equilibrium in U2 = 0.0. For this reason, it is not difficult 
to compute the Jacobian matrix for the SIR reaction term and to prove that the 
diffusion coefficients do not affect the sign of the real part of the eigenvalues. In 
other, words, the Turing space is always empty in this case (see [H]), and the 
populations tend to their constant equilibrium state for long times. 
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6.2. Example 2. An epidemic model of SARS with Patterns Formation. 

Now, we consider the following modification of the SIR model 

Mi dx\ Ami — A — a(ui,v,2, U3) — fiui, 
u 2 dx)Au 2 = a(u 1 ,u 2 ,u 3 ) - ju 2 ~ p,u 2 - H(u 2 ), 







d t ui 


- Oi 


(6.2) 


< 


d t u 2 


- a 2 






d t u 3 


- 03 



In this modified SIR model, A is the recruitment rate of the population (such as 
growth rate of average population size, a recover becomes an susceptible, immi- 
grant and so on), and TL{u 2 ) is the removal rate of infective individuals due to the 
treatment. We suppose that the treated infectives becomes recovered when they 
are treated in treatment, and 



H(v) 



r if v > 0, 

otherwise , 



where r > is constant and represents the capacity of treatment for infectives. 
We consider the same incidence term a defined in the introduction (|1.4[) . and the 
natural death rate [i is included for the recovery population. The detail about 
this epidemic model can be found in [13] 118] where the authors adopt a bilinear 
incidence rate more simple than (|1.4|) . One of the application to consider this model 
is that it supposes that the capacity for the treatment of a disease in a community 
is a constant r, in order to use the maximal treatment capacity to cure or isolate 
infectives so that the disease is eradicated [18]. It can be used for example for 
mathematical model to simulate the SARS outbreak in Beijing [19] . 

The approximation of this modified SIR model using Finite Volume is very sim- 
ilar to (|1.10p - (|1.12|) . More precisely we have 



.,,"+1 1,™ 1 \ 

, T ^ u i,K ~ u i,K ( „ \ m(a K . L ) +l 

^ m{K) At HL-US) 2. -djKJ^^ 

(6.3) K en h LeN(K) v ' 1 



i f ts\ ( A i / n+l,+ n+l.+ n+l.+ \ . n+l\ r\ 

+ m(K)\-A + <r{u 1 J K ' ,u 2K ,u 3K )+/iu 1 ^ r J=0, 



(6.4) K en h LeN(K) y ' ' 



m{K) - (7 + n)u^ H«)) = 0, 

*n{pK<L) 
d(K, L) 



n+l _ n I 
tl/ ~. U 3,K U 3,K ( ST^ n \ m \VK,L) , n +l n+1 

(6.5) K en h l^n(k) 



i(K) (lul K + U{u n 2 ^)-nu^ 



= 0. 



The results of existence of solution of the Finite Volume scheme (Proposition ^. 2p . 
nonnegativity of the scheme (Lemma I2.3|) and convergence main result (Theorem 
ll.ip can be easy generalized using straightforward calculations to this modified SIR 
model (16.211 and his Finite Volume scheme (16.3I)-(16.5I). 



24 



MOSTAFA BENDAHMANE AND M. A. SEPULVEDA C. 



6.2.1. Analysis of spatial patterns. The equilibrium state of the system (16. 2p is not 
exactly the same described in |13j because we consider here a different incidence 
term than the given by Liu which takes the reduced form a — auv. In our case, the 
two positive equilibrium points are given by E\ = (u\,V\, w\) and E% = (u2, V2, W2) 
where 



. A-r A ^y/(ra-AR -Aa) 2 -4A2R a) 

(6.6) VX,2 = "7TT> 7T- ± 



2 Rq 2 a 2aR 
A-r - RqVx.2 JVL2+r 

(6.7) Ui t2 = , Wi,2 = , 

A 4 M 
with Rq = [i + 7. The linear stability of the system (|6.2p is obtained when the real 
part of the eigenvalues of the Jacobian matrix 



J = 



u+v+w 1 (u+v+w) 7 ^ u+v+w ' (u+v+w) 2 (u+v+v^ 2 

a u a uv 



u+V+w (u+v+w)' 2 u+v+w (u+v+w) 2 ' " (tt+u+to)^ 

7 -/X 

are negative. The eigenvalues of this Jacobian matrix are given by Ai = — (j, < 
and A2 and A3 which are root of 

A 2 + ( (2/, + 7) + a^^-) X+((P + 7)M + a + = 0. 



U + V + w J \ ' u + v + w 

Real part of A2 and A3 are negative if and only if all the coefficients of this polyno- 
mial function of degree 2 are positive, that is 

„. u + v + w ( u — V U V 
(6.8) > max <^ , — 

We consider in this example the parameters A — 3, u = 0.3, a — 3.8, r = 0.5, and 
7 = 0.8 (they are the same of the example of Liu [13], except the parameter a of the 
incidence term). Replacing this parameter in (|6.6| - (|6.7p . we have the equilibrium 
state 

Ex = (7.217163781 , 0.3044098832 , 2.478426355) 
E 2 = (4.010906415,1.178843705,4.810249881). 

It is easy to verify that the equilibrium point E2 verify the linear stability condition 
(|6.8j) and E\ corresponds to a instable point. Formation of spatial patterns results 
from the diffusion-induced instability when the real part of at least one of the 
(A x \ 

eigenvalues of J = d-2 is positive [17] • If we suppose that dx = 
V A 3 J 

ds, choosing the parameters above and the the equilibrium point E2, then Turing 
stabilities appear when (5.114590203 + 9 dx) d 2 + 3.289620038 - 2.20024467 dl < 0. 

6.2.2. Simulations with local and nonlocal diffusion. Similar the above Example 1, 
we consider here two simulations, one with a constant diffusion and another one with 
nonlocal diffusion. We take the same square domain of example 1, with a Cartesian 
grid and choosing N x = N y = 300 for both simulations. The discretization in 
time is given by N t = 100 time steps for T = 2.5. That is, St = T/N t and 
m(K) — l/(N x N y ). The parameter of the SARS model are given by A = 3, 
H = 0.3, a = 3.8, r = 0.5, and 7 = 0.8. 
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FIGURE 4. Patter Formation for the SARS model (constant diffu- 
sion): (ric) Random initial condition; (a) Susceptible population; 
(b) Infected population; (c) Recovery population. 



In order to observe patterns formation, we consider the initial condition as follow: 

Ui t o(x,y) = ui +E1W1; 
U2,o(x,y) = u 2 + e 2 uj 2 ; 
u 3 ,o{x,y) = u 3 + E3UJ3; 

where u x = 4.010906415, u 2 = 1.178843705, u 3 = 4.810249881 (that is the stable 
equilibrium state E 2 ), Ej — e 2 = £3 = 0.001, and u>i G [0, 1] are random variables, 
with i = 1,2, 3. 

First we consider the model with constant diffusions, with ai — 0.1, a 2 — 0.0001, 
0,3 = 0.1 (see Figure|4j. We represent the random initial condition of the susceptible 
population on the left top of the Figure |4] (the graph of random initial condition is 
is indistinguishable to the three populations). In the same figure, we can see the 
graph of the three population at time t = 2.5. 
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Figure 5. Patter Formation for the SARS model (nonlocal dif- 
fusion): (a) Susceptible population; (b) Infected population; (c) 
Recovery population. 



On the other hand, in figure [5] we observe another simulations with a nonlocal 
diffusion. In this simulation we consider the diffusion rate coefficients given by 

M if s > M 

(().!)) 2(00 = { a if e < d l M 

yS Ui J 

z if s < e, 

with M = 10 4 , e = 10" 4 , di = d 3 = 400000 and d 2 = 400. The choice of d, is 
in order to take diffusion coefficients close to the values of the constant coefficient 
diffusion of the first simulation (Figured]) a t time t = 0. In this case of nonlocal 
diffusion, the coefficient diffusions are not constant in time, obtaining obviously 
different results than the first simulation of this example 2. Moreover, we remark 
very different behaviours between both simulations (compare Figure |4] and Figure 
EJ. Outside the nonlocal diffusion, we consider the same parameters for both sim- 
ulations in this example 2. Finally, we observe in Figure El the evolution in time 
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Figure 6. Example 1: evolution in time of the diffusion terms 
a i (^J u i ' w ith i = 1,2,3; Left: Susceptible and recovery pop- 
ulations; Right: Infected population. 



of the diffusion at 
diffusion. 



Ui dx I , with i = 1,2,3 for the simulation with the nonlocal 



Appendix A. Uniqueness of weak solutions. 

In this appendix we prove uniqueness of weak solutions to our systems by using 
duality technique (see e.g. [H]), thereby completing the well-posedness analysis. 

First, we consider (u%, u 2 ,u 3 ) and (i>i,i>2,V3) two solutions of the system (|1.1|) - 
(|I.3[) . We set U = Ui — w, for i = 1,2, 3, then U satisfies 

dtUi—(ai( Ui dx^j Aui — a>i ( / vidx^jAvi 

= -(a(u 1 ,u 2l u 3 ) - <j(vi,V2,v 3 )) - fiUi, 

dtU 2 -(a 2 (^J u 2 dx^JAu 2 -a 2 ^J v 2 dx^JAv 2 

= (<t(ui,U2,u 3 ) - a(vi,v 2 ,v 3 )) - + fi)U 2 , 

d t U 3 -(a 3 (^J u 3 dx^JAu 3 -a 3 (^J v 3 dx^j Av 3 ^j = jU 2 , 

Vui ■ r\ — Vvi ■ r\ — on Sx, i = 1, 2, 3, 
Ui(x,0) = for x e Q, i= 1,2,3. 

Now, we define the function tpi solution of the variational problem 



(A.l) 



(A.2) 



Vipi(t,-) -V(f)dx = / Ui(t,-)(j>dx, 
iq Jn 
for all 4> S H (fl), such that J n <j)dx 



$ = 1,2,3, 







for a.e. t 6 (0,T). Since Ui and Vi are in L 2 (Qt), then we get from the theory 
of linear elliptic equations, the existence, uniqueness and regularity of solution tfi 
satisfying 

ifi £ C([0, T]; JJ 1 (f2)) with / ^-(i, •) dx = 0, for i = 1,2,3. 
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Note that from the boundary condition of ipi in (IA.2j) and Ui(Q, •) = we deduce 
that 



(A.3) Vv3i(0, •) = in L 2 (ft) for i = 1, 2, 3. 



Multiplying the first, second and third equations in (jA.ll) by i()i,ip2,ip3 E L 2 (0, T; i? 1 (Sl)), 
respectively, and integrating over Q t := (0,i) x fi, we get 



(A.4) 

/ {d s Ui, ipi) ds + ^ai ( I u\ dx^j Vu\ — a i(l Vi dx\ Vvi^j -Vipi dx ds 
^—■y J o <j <j q± «/ o «/ ri 

+ (a 2 (^J^ u 2 dx S jVu 2 - a 2 (^J v 2 dx^j\7 v 2 ^ ip 2 dx ds 
a 3^J u 3 dx^J V«3 — 03 U3 da;^ V«3^ -VV>3 da; 
(cr(ui,U2,U3) - <t{vi,v 2 ,v 3 )) - LiUi^ipi dxds, 

2^3 dx ds. 



j -02 dx ds + 


S Sq 







Since ifi £ £ 2 (0, T; we can take ipi = ipi in (IA.4|) and we obtain from (I A.2|) 

and ((A3)) 



2 / (d s U i: (pi)ds = -2 / (d s A(pi,ifi}ds 
Jo Jo 

(A.5) = / |V^(t,x)| 2 dx- / |V<^(0,a;)| 2 da; 



|Viy9i(t, a;)| 2 dx. 



£2 



From definition of a we obtain easily 



(A. 6) \a(ui,u 2 ,u 3 ) - a(vi,V2,v 3 )\ < |«i - v t \ , 



for some constant C > 0. 



FINITE VOLUME SCHEMES FOR NONLOCAL REACTION-DIFFUSION SYSTEMS 29 



Using (|A.2|) . (IA.6j) . Holder's, Young's, Sobolev poincare's inequalities yields from 
ljA.4jl with ip t = ifi 

(A.7) 

/ {d s Uj, ipi) ds = — // Off / Uidx)UiAipidxds 
i=1 Jo i=\JJQt y ->& ' 

— // f a «(/ u idx\—cii(^ Vidx\Vvi ■ Wtpidx ds 



^(cr(ui,M 2 ,u 3 ) - v(vi,v 2 ,v 3 )) ~ nU^jtpi dxds 
+ II (^((t(u 1 ,u2,u 3 ) - a(vi,V2,va)) - (j + (j,)U2^ <P2 dx ds 
jU2<f3 dx ds 



-VM, // \Ui\ 2 dxds 
t=i JJ Qt 

3 ft/f, r r ^ /.t 

+ JJ Q \ u i\ 2 dxds + C iYlJ o W Wv i\\l 2 (n) W Vl Pi\\ 2 L2 {n) ds 



3 

+ ^ 



C 2 J o WVipxWl^ ds + Yj^ JJ^ \U t \ 2 dxds + C 3 J^ || 
^ |[/ 2 | 2 dxds + C 4 y || V<£ 3 ll^n) ds 



M2 

12 



<(d + C 2 + C 3 + C 4 ) ^ A||V«i||i a(n) + 1) ||V^||* 2(n) ds, 

for some constants C\, C2, C3, C4 > 0. 

Using Vfj <E L 2 (Qt) for i = 1, 2, 3 and Gronwall's lemma to conclude from (|A.7|) 

= 0, i= 1,2,3, 

almost everywhere in Qt, ensuring the uniqueness of weak solutions. 
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